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Abstract 

First we review the derivation of the relativistic flux tube model for a quark- 
antiquark system from Wilson area law as we have given in a preceding paper. 
Then we extend the method to the three-quark case and obtain a Lagrangian cor¬ 
responding to a star flux tube configuration. 

A Hamiltonian can be explicitly constructed as an expansion in 1/m^ or in the 
string tension a. In the first case it reproduces the Wilson loop three-quark semirel- 
ativistic potential; in the second one, very complicated in general, but it reproduces 
known string models for slowly rotating quarks. 
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1 Introduction 


In a recent paper we have shown that, for the quark-antiquark system, ne¬ 
glecting spin-dependent terms, it is possible to derive rigorously the so-called 
relativistic flux tube model from the Wilson area law in QCD. 

To achieve the result, the essential ingredients were the path-integral represen¬ 
tation of the Pauli-type propagator we used in the derivation of the semirela- 
tivistic potential, the replacement of the 1/m? expansion of the purely kinetic 
terms by their closed form and the explicit integration over the momenta of 
the functional integral. 

In this letter we want to show how a similar result can be worked out for 
the three-quark system starting from the corresponding representation for 
the three-quark propagator we have considered in ref. [1]. The resulting La- 
grangian turns out to correspond to a star flux tube conhguration. As in the 
quark-antiquark case, the Hamiltonian can be expressed only in the form of an 
expansion in 1/m^ or in the string tension a. In the first situation we obtain 
the Wilson loop semirelativistic potential, but for the spin-dependent part. 
In the second one we end up with a very complicated result which, however, 
reduces to the three-quark string model already considered in the literature 
for low angular momentum. 

In Section 2 we review the quark-antiquark case to clarify the procedure and 
to establish notation; in Sec. 3 we derive the new Lagrangian and in Sec. 4 we 
obtain the Hamiltonian of the resulting model. 


2 Quark-antiquark flux tube model 


In QCD the quark-antiquark gauge invariant propagator can be written 


G{xi,X2,yi,y2) = 

= ^^^{U{x2,Xi)Sl{xi,yi\A)U{yi,y2)C~^Sl{y2,X2\A)C) . ( 1 ) 

Here c denotes the charge-conjugate helds, C is the charge-conjugation matrix, 
U the path-ordered gauge string 


/ 

U{b, a) = P exp ig j dx^ 


[X] 


( 2 ) 
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(the integration path being the straight line joining a to b), Sf and the 
quark propagators in an external gauge held A'^. Furthermore, in principle the 
angular brackets should be dehned as 


JVAMf{A)e^s[A] ’ 


( 3 ) 


5[A] being the pure gauge held action and Mf{A) the determinant resulting 
from the explicit integration on the fermionic helds. In practice in this paper 
we take Mf{A) = 1 (quenched approximation). 

From (1), by a Foldy-type transformation and with the other elaborations 
described in ref. [1], after setting Hi = 1/2 = h, Xj = X 2 = U, we obtain the 
following Pauli-type two-particle propagator 


7t'(xi,X2;yi,y2;tf - U) = 


XI 


X2 


J Pzi'Dp2 J 'Dz2pp2exp<i 

yi y 2 ^ 




dtJ2 (pj ■ % - + • • • 


i=i 


2mj 8m| 


—i In Wgg + spin dependent terms 


( 4 ) 


Wqq being the so-called Wilson loop integral 



Here the integration loop F is assumed to be made by the quark world line 
Fi, the antiquark world line F 2 (described in the reverse direction) and the 
two straight lines which connect xi to X 2 , y 2 to yi and close the contour (see 
Fig. 1). As usual A^{x) = 4AaA“(a;), P prescribes the ordering of the colour 
matrices (from right to left) according to the direction hxed on the loop, and 
the “functional measures” are given by 

3N . 3N 

( TYh \ 2 / XS \ ^ 

j dzi...dzN-i, ^ ( 2^ j dpi... dpw-idpAT, (6) 


where e = t/N, the limit N ^ 00 is understood and the end points x and y 
stand for the conditions zq = y, zn = x. 

After neglecting the spin-dependent terms and reintegrating the kinematical 
terms in their closed form, (4) becomes 
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Xl X2 

7^(xi,X2;yi,y2;tf “ ^i) = /'^zi'^P2 J 'Dz2'Dp2 x 

yi y 2 

if 2 


exp < * 


f dt^ (pj ■ Zj - ^Jm] + p|) - i In VP, 
^ i=i 


> . 


( 7 ) 


As usual we assume that i In Wqq can be written as the sum of a perturbative 
term (at the lowest order) and an area and a perimeter term 

i In J dx^ J dx^ iD^^{xi - X 2 ) + aS'^m + CP , (8) 

Ti r2 


S'min being the minimal area enclosed by F, P its perimeter and the usual 
gluon propagator. Furthermore, to evaluate S'min, we adopt the straight line 
approximation, consisting of replacing the minimal area with the area spanned 
by equal time straight lines joining the quark and the antiquark 


S'min = dt ds\zi - Z 2 I 1 - (SZIT + (1 “ s)z 2 t) 


( 9 ) 


having set ijtp = (transversal velocity), /r and r = zi — 

Z 2 . We also use the Coulomb gauge and adopt the instantaneous approximation 
on the perturbative term. So we can write 


3 r 


1 _ ^ 


+ 


2^1/2 


Par ds 1 - (sziT + (1 - s)z 2 t) + C 



( 10 ) 



Notice that in the hrst term in (8) and in the last in (10), we have neglected 
the contributions coming from the end straight lines {t = tf and t = f,) having 
in mind the limit tf — t, ^ oo. Furthermore, (9) not being Lorentz invariant, 
our assumptions are understood to be made in the centre-of-mass frame. 


Then we replace (10) in (7), expand the kinetic terms around pj = rnjZjj 


p ■ i — J + p2 = —mVl — _ z^)2((5'^*^ — z^z^) 

^ 2m 

mz^ \ f k \ 

7r^ J V ~ 7r^ j ^ ' 


(11) 
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and perform explicitly the integration over the momenta neglecting the suc¬ 
cessive terms in (11) (Gaussian approximation). We obtain 


it'(xi,X2;yi,y2;tf - h) = 

Xl X2 

J VziA[zi] J 'Dz2A[z2]exp 

yi y 2 


dtL{Zi,Z2,Zi,Z2) 


ti 


( 12 ) 


with 


L = 


i=i 

1 

ar J ds 
0 




1- - ((5“ + f^r^)ii*i2^ 


+ 


1 - (SZIT + (1 - s)z2t)' 


nl/2 




cj:j dt^i 

1 = 1*; 


J ’ 


(13) 


and 


A[z] 



_(l_z2)l/2(^hfc_^/i^fc^ 

_m 



= 11 ( 1 - 2 ^) ( 14 ) 

t 


Of course the factor A);^] has to be considered part of the relativistic “func¬ 
tional measure” in the conhguration space. 

Notice that if we introduce explicitly the factor l/% in front of the exponen¬ 
tial in (7), the result expressed by (12) and (13) becomes exact in the limit 
h —* 0. This implies that (13) provides already the exact classical Lagrangian 
of the system, while to go beyond the Gaussian approximation would amount 
to modifying the expression of A(z) alone. As a matter of fact, the addi¬ 
tional terms would be highly singular in the time lattice spacing e and would 
not match in reconstructing a time integral in the exponent (the fact can be 
checked explicitly performing the exact integrals in the momenta [2]). 

The Lagrangian dehned by (13) is the Lagrangian of the relativistic flux tube 
model [3,4] with the addition of a Goulombic term. The perimeter term can 
be absorbed in a redehnition of the quark masses {rrij —* rrij + C) and from 
here on we ignore it. 

From (13) we can introduce the canonical momenta 


Pi 


dL 
dir I 


miVi 

\/l - vf 


-|- ar 


ds 


0 


-SVt 
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P2 = 


dL 

dwo 


m2V2 , / . (1 - '5)vt 

+ ar I as — , 

1-v? 


1-vi 


0 


(15) 


while the total linear momentum and the Hamiltonian turn out respectively 
to be 


P = Pi + P2 = 


miVi 7712^2 , f , Vt 

H- 1 +ar ds- 


vf 


vi 


V? 


(16) 


H = Y.Vr^j-L = Y. 1^ 2 

j=l ^|l - 


i=i 


worlds- ^ 


1-v? 


(17) 


Here we have set Vj = Zj and by Vt(s) = szit + (1 — s)z 2 t we have denoted 
the velocity of the flux tube segment specihed by s and s + ds. 

Eqs. (15) cannot be inverted in closed form, but it can be by an expansion 
in 1/m^ or in a. In the hrst case we reobtain the semirelativistic Wilson-loop 
potential [1] (apart from the spin-dependent terms), in the second one we hnd 
[4] at the hrst order in a 


H = ^ml + + 


ar 


^Jmi + + ^Jm2 + 


X 


\Jml + q^^Jml + / 


arcsm 




mf -|- q^ 


+ arcsm 


gx 


+ g^ 


(18) 


where relevant simplihcations have been obtained setting explicitly pi = 
—P 2 = q (centre-of-mass frame) and where we have dehned Qr = (q ■ h)f . 

Notice that if we use (18) in the phase space path integral and try to go back 
to Eq. (12) by integrating again over the momenta, we would hnd a diherent 
A(z) due to the occurrence in (18) of a momentum-dependent interaction part. 
This means that to be consistent we have actually to introduce an appropriate 
normalization factor in front of the expression of Wqq as would be given by 
(8) (see [7]). 
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3 Three-quark flux tube Lagrangian 


The three-quark gauge invariant propagator is 


G{xi,X2,X3,yi,y2,y3) = 

_ 1 

(0|T {xm, {xm, X2)U°‘^''^ {xm, a;i)^/’3c3 (a^3)^2c2 {x2)i^ici (a^i) 

1 


21 ^aia2a3^bib2bi 

U{xm,xi)SI {xi,yi\A)U{yi,yM) 
U{xm, X2)S2 {x 2, y2\A)U{y2, yu) 
u{ xm: X 3 )S^{xs, y 3 \A)U{y 3 , yu) 


aibi 


02^2 


azbs 


(19) 


where we assume x^ = x^ = x^ = = tf, yi=y 2 =y 3 = y^M — and 

jM are points chosen inside the triangles (xi, X2, X3) and (yi, y2, ys) in such a 
way that Z]j=i \^j ~ ^m\ and Z]j=i |yj ~ YmI are minima and a*, bi are colour 
indices. 


Again, neglecting spin-independent terms, we can write 


A'(xi,X2,X3;yi,y2,y3;tf - hi XM,yM) = 

XI X2 X 3 

= J VziVpi J VZ 2 VP 2 J VZ 3 VP 3 

yi y 2 ya 


exp 


y 3 

dtY, (Pi ■ % - 



i In W 3 q 


•> 


( 20 ) 


I L‘i JJ 

which corresponds to (7). W 3 q is the “Wilson loop integral” for three quarks 
defined by 


1 

/ 


( f . 

ff"3a = -r' 

iq 3 

y ^010203 ^616263 

P exp 

ig J dx^^A^^{x) 

\ Ti / 


" 

( f 

0262 

“ 

( f \ 

P exp 

ig / dx>^^A^^{x) 


P exp 

ig / dx^^^Af,^{x) 


\ r 2 / . 



\ T 3 / 


( 21 ) 


as 63 


1 
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Fj being made by the world line Fj of the quark j plus the straight lines 
connecting i/m to i/j and Xj to xm (see Fig. 2). As in (8) we shall set 

i In Wsq = j - Xj) + aS'min + CP , (22) 

^<Jr* r 

i-1 1 j 

where the hrst term is again the lowest-order perturbative contribution; P is 
the total length of Fi, r 2 and Fs, S'mm is the area of the minimal three-sheet 
surface enclosed by Fi, r 2 , Fs and a mid-point world-line Tm {zm = ZM(t)) 
joining yM and Xm ] the minimum is performed for hxed quark world line, but 
varying Tm- 

Again adopting instantaneous and straight line approximations we obtain 


i In Wsq = 




-2 
3 Tqi 




h \ 'h 'k 


+ 


3 


i=i 


1 



0 





(23) 


Here = Zj — Zj, yj = — zm and zm has to be chosen in such a way that 

rj is minimum at any given time, Zjrj-. = {5^^ course, ZM(t) 

has to be coplanar with zi(t), Z 2 (t) and z^it) and then two different types of 
conhgurations are possible: 


I) if in the triangle (zi(t), Z 2 (t), z^it)) no angle exceeds 120°, then ZM(t) has 
to be such that ri(t), r 2 (t), r 3 (t) make angles of 120° with each other; 

II) if the angle in zi{t) reaches 120°, then zm = z;- 


Notice that correspondingly for zm we have 


Zm — 


z; 


I type conhguration 

II type conhguration 


(24) 


TZ being the matrix with elements — f'jf^)/rj. 

Replacing as before (23) in (20) and integrating on the momenta in the Gaus¬ 
sian approximation we can write 
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J=1 


Notice that again the perimeter term can be absorbed in a redehnition of the 
quark masses. 


4 Three-quark relativistic flux tube Hamiltonian 

For simplicity, in this section we shall omit the perturbative part of the 
Lagrangian, keeping in mind that this can be simply added to the pure 
flux tube result. To simplify the notations let us set v* = zp vm = zm, 
= 'SjVjj'. + (1 — Let us also introduce the matrices 7^ with 

elements 7*^^ = <5^^ — then we have obviously 

n = (27) 

j=l ' 3 

and 


ViTi = 71 Vi, wm = 'R > 

i=i 

^ 1 

v*(si) = 71(siVi + (1 - Si)vM) = SiTiVi + (1 - Si) ^ — 717^ v^ . (28) 

j=i 
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So we can write the flux tube Lagrangian as 


L = 


i=l i=l Q 


( 29 ) 


and, we can obtain the canonical conjugate momenta 


rriiWi 


Pi = 




ds 


{sjdij + (1 - Sj)TZ ^Ti/r-i) 


i=i 


1 - V 


t2 


(30) 


while the total linear momentum and the Hamiltonian turn out, respectively, 
to be: 




i=l 


=1 a/ 1 - v/ 


2=1 


1 - V 


t2 


(31) 


^ rrii 


H = J2Pi-^i - L = 


3 


i=l 


i=l 


+ ay r-i dsi 


1 Vj l—\ g 


1 - V, 


t2 


(32) 


Again Eq. (30) can be solved with respect to Vj (j = 1,2,3) only as an 
expansion in l/m^ or in the string tension a. Consequently even H turns out 
to be expressed in such a form. 

Using an expansion in l/m| we obtain 


3 1 3 

i=l ^i=l 



d2 1 3 
Pj _ V~' Pi 

s 2-^ 


m. 


i=l 


rrii 


, _2 , PiTi ■ VMTi 

+ ^MTi ' _ 

mi 


+ 0 


Vi=Pi/™i 



(33) 


with 


VMTi 


'fi='Pilrr>.i 






PjT, 

mj 


(34) 


Eq. (33) coincides with the static plus the velocity-dependent part of the 
three-quark Wilson loop Hamiltonian [1]. 
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Considering an expansion in a in (30) we obtain the explicit relativistic Hamil¬ 
tonian 


^ = E 

2=1 


pf + 771^ + a 


3 

2=1 


i j dsi \J~^ 


,t2 


+ 0(a‘^) ,(35) 


Vi=Pi/^/pf+mf 


obviously now 


* IVi = 


= Si 


P*T, 




mf -h p: 


+ (1 — Si)TZ 


3 1 

-1 _T. 


PjT, 


i=i 


J 


ml -h 


(36) 


Notice that like Eq. (9), even Eq. (23) and so even Eqs. (26), (29), (32) are 
supposed to hold in the centre-of-mass frame P = pi -|- p 2 -|- Ps = 0 and 
Pi,P 2 ,P 3 are not independent. A possible choice of independent variables is 
given by the internal Jacobi momenta, in terms of which we can write 


Pi = q, 


P2 = 


m2 


m2 + ms 


q + k , Ps = 


ms 


m2 + ms 


q-k. (37) 


Eq. (37) should be used in (36), (35). The hnal result is very complicated. Eor 
low angular momenta, however, nP should be negligible and Eq. (35) reduces 
to the three-quark starlike string model 


^ = E 


2=1 


+ Pi + a 

2=1 


(38) 


already considered in the literature [5,6]. 

Performing explicitly the s* integration in (35) we obtain 


H = J2yPi oA/i -+ 

i=l i=l I ^ 


VjUTi ■ (VjTi — ^MTi) 


2(ViTi - 
1 


1 - V?T,. 


1 — 


2|VjTi — ViWTil 


X arcsm 


(1 + 


(VjTi ■ (ViTi - VmtJ)" 


(ViTi - VmtJ^ 
ViTi ■ (ViTi - VmtJ 


arcsm 


^(1 - - VmtJ^ + (ViT, ■ (ViT, “ VmtJ)^ 

VMTi ■ (VjTi — VmtJ 


\/(l - ^MTi)(yiTi - VmtJ^ + (ViT, ' (VjTi “ VmtJ)^ 


>(39) 
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which, however, does not essentially simplify after the substitution (37), con¬ 
trary to what happens in the quark-antiquark case (18). 


5 Conclusion 


In conclusion we have shown that even for the three-quark system a flux-tube 
like relativistic Lagrangian can be obtained in the Wilson-loop framework, if 
one neglects spin. The method does not seem appropriate for a direct intro¬ 
duction of spin since it uses a Foldy-Wouthuysen type transformation at an 
early stage and so the spin dependence in the starting path-integral occurs as 
an expansion in l/w?. 

As a matter of fact, in the quark-antiquark case the relativistic flux-tube model 
is strictly related to the instantaneous approximation of the corresponding 
Bethe-Salpeter equation, and spin can be taken into account through such an 
equation [8]. In principle the generalization of such an equation should be the 
convenient formalism even for the three-quark bound state. 

Notice finally that to give the quantum Hamiltonian operator corresponding 
to Eqs. (35) and (39), attention has to be paid to questions of ordering. If the 
discrete form of the quantity W^q is written as [1] 


21 ( ^aia2a3^f)lfe2b3 


(40) 


n 

t=i 


P exp ^ (4n - 4n-l) ( 


n T n—1 


the correct ordering would be the Weyl’s ordering. 
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Fig. 1. Quark-antiquark Wilson loop 
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Fig. 2. Three-quark Wilson loop 
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